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The dynamics of stars around spiral arms 
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ABSTRACT 

Spiral density wave theory attempts to describe the spiral pattern in spiral galaxies in 
terms of a long-lived wave structure with a constant pattern speed in order to avoid the 
winding dilemma. The pattern is consequently a rigidly rotating, long-lived feature. 
We run N-body simulations of a giant disc galaxy consisting of a pure stellar disc and 
a static dark matter halo, and find that the spiral arms are transient features whose 
pattern speeds decrease with radius, in such a way that the pattern speed is almost 
equal to the rotation curve of the galaxy. We trace particle motion around the spiral 
arms. We show that particles from behind and in front of the spiral arm are drawn 
towards and join the arm. Particles move along the arm in the radial direction and 
we find a clear trend that they migrate toward the outer (inner) radii on the trailing 
(leading) side of the arm. Our simulations demonstrate that because the spiral arm 
feature is co-rotating, the particles continue to be accelerated (decelerated) by the 
spiral arm for long periods, which leads to strong and efficient migration, at all radii 
in the disc. 

Key words: galaxies: evolution - galaxies: kinematics and dynamics - galaxies: spiral 
- galaxies: structure 



1 INTRODUCTION 

Spiral density wave theory, fi rst postula t ed by iLindbladI 
(|l96d l and carried forward by iLin Shul (|l964l l , has been 
the most widely accepted theory of the spiral arm for al- 
most 50 years. This theory considers the spiral arm to 
be the crest of a stellar density wave that rotates around 
the centre of the Galaxy at an angular pattern speed, lip, 
^ that does not vary with galacto-centric radius. Spiral den- 
sity wave theory naturally explains the long-lived nature of 
the spir al arm. However, short-lived spi r al arms were pro- 
posed bv lColdreich Lvnden-Belll (|l965l l. I Julian k Toomr3 
(|l966l l and lToomrd (|l98lh . These studies argue that spiral 
structure can grow from local instabilities and shearing ro- 
tation predicted from a flat rotation curve of disc stars that 
help to develop the spiral features (|Athanassoulalll984l V Al- 
though these spiral arms are short-lived features that con- 
tinuously disappear, they are regenerated by local growth 
of instability. In these transient spiral theories, the spiral 
structure itself can be self-excited in small oyer- d ense re- 
gions as a consequence of shot noise (|Toomrelll99Qi ). If the 
over-density increases above a few percent, the growth be- 
comes very non- linear and results in rapid growth of the 
spira l arm density (e .^. Toomre 198ll: lDonner T homassonI 
1 19941 : lBotte"^ l20Q3l: Sellwood 2010). Non-linear growth of 
transient rigid waves is also demonstrated in the context of a 
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shearing sheet (jPuchs et al.ll2005l l. While shearing is impor- 
tant for non-linear growth of spiral arms, it is also reported 
to be the main fac tor in the a ppearance of spurs and spiral 
arm bifurcation (|Dobbs Bonnellll2008l ). 

More recently, the development of high-resolution nu- 
merical N-body/SPH simulations has played an important 
role in the study of spiral structure formation and evolution. 
Many numerical simulations including self-gravity and with 
no addition or forcing of spiral potential show short-lived 
transient spiral arms, which come and go throughout the 
evolution suc h that the spiral pattern o f the disc is never 



removed (e.^. ISellwood Sz Carlberd 1983: iBaba et al 



2OO9I : 



"Quill en et al.1 l201ll:lF\liii et al.1 I2OIII : IWada et al.1 l201lh . 

Sellwo od Carlberd ([Tgsi r note that the spiral pattern in 
N-body simulations generally fades over time because the 
spiral arm structure heats the disc kinematically and causes 
Toomre's Q parameter to rise, hence the disc becomes too 
stable against the development of non-axisymmetric struc- 
ture. They suggest that continuous addition of a kinemati- 
cally cold popul ation of stars is crucial t o ma intain the spi- 
ral arm pattern. ICarlberg FreedmanI (|l985l ) demonstrate 
this with gas clouds in circular orbit that are allowed to 
become future sites of star formation. Recently iFuiii et ahl 
(2011) suggest that the rapid disappearance of the pattern 
probably results from a low number of particles in these 
simulations, which leads to higher Poisson noise and more 
rapid growth of strong spirals, which then heat the disk and 
erase the spiral pattern. They demonstrated that the spiral 
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pattern can be sustained without dissipation or continuous 
addition of cold components if a sufficiently high number of 
particles (> 10^) is used to prevent a large Poisson noise 
and therefore fast disc heating. Still, there are no numerical 
simulations that reproduce long-lived s piral arms like those 
predicted by spiral density wave theory (|Sellwoodl[20lTl l. and 
this is the case fo r both spiral and barred-spiral galaxies (e.g. 
iBaba et al.|[2009l : iQuillen et aLll201ll : IWada et al.l l201lh . 

In spiral density wave theory, the angular pattern speed 
of the spiral arm is constant as a function of radius. The 
flat rotation curve of spiral galaxies implies the decreasing 
angular velocity of stars. Therefore, there is one radius 
where the stars rotate at the same velocity as the arm. This 
is called the co-rotation radius, inside which stars rotate 
faster than the spiral arm and outside which they rotate 
more slowly. The co-rotation radius marks a point of reso- 
nance, where stars very close to this radius have been shown 
to undergo large changes in their angular momentum and 
migrate radia lly to another part of the disc . This was first 
proposed by iLvnden-Bell Kalnaid (|l972 ). who showed 
that radial migration around co-rotation ca n occur without 
increasing the random motions of stars. Sel lwood BinnevI 
(I2OO2I ) extended this theory, and showed that transient 
arms provide a good mechanism for mixing, where each 
transient arm possesses a unique co-rotation radius around 
which nearby stars can migrate. The main conclusions of 
this paper were that only around the radius of co-rotation 
is there large angular momentum transport of individual 
stars without heating. There have been many studies on the 
effect of co-rotation resonances on str ucture, radial heatin g 
and r aixing in spiral g a laxies (e.^. [ Carlberg SellwoodI 
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T997I: ISellwood BinnevI I2OO2I: 



QuillenI 120061: iRoskar et al.ll2008la: iRoskar et al 



anchez-BlazQu ez et al 
and in barred-spiral ga laxies where 



iMinchev et al 
bar 



nance is also include d (e.g. IMinchev FamaevI I2OIOI : 

iManos Athanassoulall201lh . 

Very recently, there have been observational studies 
with evidence of a pattern speed that varies with radius, 
which is different from the constant pattern speed of spiral 
density wave theory. The evidence comes from a range of 
techniques and sources that include observations of no off- 
set between different s tar formation tracers in nearby galax- 
ies (Fovle e t al.ll201ll ). and radially decreasing spiral arm 
patt ern speeds deduc ed by the Tremaine- Weinberg method 
(e.g. iMerrifield et al.1 2006; M eidt et al.1 12OO8I : iMeidt et al.1 
l2009l :ISDeights WestpfahLi201ll ). which inevitably leads to 
short-lived arms as a con sequence. Numerical simulations of 
an interacting galaxy in Dobbs et al.l (|2010) also shows a 
pattern speed that decreases as a function of rad ius. During 
the preparation of this work, IWada et al.l (|201ll ) performed 
simulations of an isolated galaxy, and found a pattern speed 
that traces very well the rotation curve, which indicates that 
there is co-rotation of the spiral pattern with the stars at all 
radii. 

In this study, we analyse the pattern speed of the spi- 
ral arm feature and the dynamics of particles around the 
arm using our N-body simulations of a pure stellar disc. We 
determine the nature of the spiral arm pattern speed, and 
monitor how stellar motion is affected by the arm. We also 
find that the spiral arms in the simulations are transient 
and that the pattern speed is always similar to the circular 
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Figure 1. The initial circular velocity (dot-dashed red line) and 
the circular velocity at t= 1.77 Gyr (solid black line). 



velocity. We offer a qualitative insight into how phenomena 
such as radial migration might occur with the transient co- 
rotating spiral arms, and how the co-rotating spiral arm may 
evolve. 

In Section 2 we explain how we set up the model and the 
initial parameters that we choose. In Section 3 we present 
the results of our analysis, compare them with previous stud- 
ies and discuss their implications. In Section 4 we summarise 
the significance of the results and remark upon the value of 
the simulations and future work. 



2 METHOD AND MODEL SETUP 

Our sim ulations are performed with a Tree N-body code, 
GCD+ (jKawata GibsonI [20031 ) . We set up a disc galaxy 
that consists of a dark matter halo and a pure stellar disc 
with no bulge, and which is similar (slightly smaller) in size 
to the Milky Way. We describe the stellar disc component 
as a collisionless N-body system of particles, and adopt a 
static dark matter halo potential. 

The dark matt er halo density profile follows that of 
iNavarro et al 

where Pc is the characteristic density described by 
Navarro et al. the concentration parameter, c — 

and X = "^^^ scale length is rg, and r2oo is the 

radius inside which the mean density of the Dark Matter 
sphere is equal to 200y9crit: 

M200 



r2oo = 1.63xl0"^( 







*[7^^]^(l+^o)-^/l-^kpc.(2) 



We assume M200 1.7 x 10^^ M©, c = 15, Ho 0.266, 
2:0 = and Ho — 71 km s~^ Mpc~^. We do not apply 
adiabatic contraction for the dark matter halo, for simplicity. 

The stellar disc is assumed to follow an exponential sur- 
face density profile: 



Md 1 2/ ^ N f R \ 



(3) 
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Figure 2. The ratio of velocity dispersions in the radial and z 
direction, fji = (Jr/(Jz, at t=0 (dot-dashed red line) and t = 1.77 
Gyr (solid black line), plotted as a function of radius. 
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Figure 3. Toomre's instability parameter, Q, computed at t = 
(dot-dashed red line) and t = 1.77 Gyr (solid black line), as a 
function of radius. 



We apply the mass of the disc, Md = 3 x 10^° M©, the scale 
length, Rd = 3.0 kpc and the scale height Zd = 0.35 kpc, 
which is constant over the disc. These parameters lead to a 
stellar surface density of E = 37 Mq pc~^ at 8 kpc, which is 
similar to S = 35.5 Mc?) pc~^ obtained for disc stars in the 
solar neie: hbourhood (|Flvnn et al.ll2006l ). 

The initial rotation curve is shown in Fig. [T] The rota- 
tion speed at 8 kpc is 210 km s~^, which is slightly lower 
than the ge nerally accepted value of 220 km s~^ or even 
higher (e.g. iMcMillanI l201lh for the Milky Way. Although 
we constructed a galaxy similar in size to the Milky Way, 
it is not the intention of this study to reproduce the spiral 
arms of the Milky Way. We also deliberately choose a disc 
galaxy model that does not develop an obvious bar structure 
in order to avoid complexity resulting from the bar poten- 
tial, hence we can focus more on the pure effect of spiral arm 
development. iFuiii et al 1 (|201lh show that if more than one 



million particles are used to describe the disc component, 
artificial heating that suppresses the spiral arm formation is 
not significant. We use 3 x 10^ particles for the disc com- 
ponent, therefore the mass of each particle is 10^ M0. We 
adopt a fixed softenin g length of 340 pc , with the spline 
softening suggested bv lPrice Monagha 3(|2Q03). 

The velocity dispersion for each t hree dimensional posi - 
tion of the disc is computed following "Springel et alj (|2005t ) 
to construct the almost equilibrium condition. One free pa- 
rameter in this method is the ratio of the radial velocity 
dispersion to the vertical velocity dispersion, /h, which re- 
lates as fn = (JrIoz- We choose /k = in the simulation 
shown. This is slightly lower than . fa ^ 2 that is the ob- 
served ratio for the Milky Way (e.g. iHolmberg et al.l [20091 : 
iBinnevllioTol l. However, as mentioned above, we do not aim 
to create a Milky Way- like galaxy. 



3 RESULTS AND DISCUSSION 

The simulation set up in Section 2 is evolved for 2 Gyr. As 
with the previous studies described in Section 1, we also find 
that the disc develops transient and recurrent spiral arms. 
In this paper, we focus on one arm identified around i — 1.77 
Gyr. The circular velocity at t = 1.77 Gyr is shown in Fig.[Tl 
which is not significantly different from the initial circular 
velocity. Fig. [2] shows fn as a function of radius at t = 
and t = 1.77 Gyr. The value drops slightly with time as the 
disc heats up slightly during evolution. This is quantified 
in Fig. O which shows a slight growth of Toomre's instabil- 
ity parameter, Q, at the same time step. However, there is 
not significant he ating, which helps to maintain spiral arms 
(|Fuiii et al.ll201lh . 

Although we mainly discuss the analysis around i — 
1.77 Gyr in this paper, we also applied similar analyses to 
other spiral arms that developed at different times in this 
simulation as well as spiral arms in other simulations with 
different initial configurations of the disc and dark matter 
halo. We find that all the spiral arms we analysed show very 
similar results to those shown in this section. We also find 
similar results in several (barred and non-barred) simula- 
tions that take gas and star formation into account. These 
will be described in a forthcoming paper. 

First, we present the analysis and results of the pattern 
speeds of the chosen spiral arm. Then we examine the mo- 
tion of selected particles around the arm, and present and 
discuss an analysis of their angular momentum and energy 
evolution. 



3.1 Pattern Speed and Radial Migration 

The middle row of panels in Fig. |4] shows three consecutive 
snapshots of the model galaxy, where the white lines depict 
the peak density of the chosen spiral arm at each radius. This 
is calculated by first creating a smoothed normalised den- 
sity distribution in the R-Q plane, as in the left and middle 
columns of Fig. [8] At each radial bin, we locate the azimuth 
of maximum density. In this way, we traced the arm in the 
radial range between 4.5 and 8.5 kpc, and so we focus our 
analysis on this radial range. The pattern speed, Qp, is then 
easily found from the azimuthal offset of the peak densities 
between the time steps. 
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t=1.71(Gyr) 



t=1.73(Gyr) 






t=1.75(Gyr) 



t=1.77(Gyr) 



t=1.79(Gyr) 






t=1.81(Gyr) 



t=1.83(Gyr) 



t=1.85(Gyr) 



Figure 4. Snapshots of the disc from face on view. The white hne marks the position of the peak density of the fully formed spiral 
arm at the time indicated (see text for how this is determined). In the middle row, the white crosses at 5 kpc and 7 kpc radius indicate 
position of the peak density line at those radii. In the top and bottom rows, the crosses have been rotated from their positions in the 
middle panel (t= 1.77 Gyr) by an angle derived from the angular rotation speed in Fig. [5] for those radii. The purpose of the crosses is 
to guide the eye to the spiral arm from formation around t= 1.73 Gyr, to its apparent destruction around t= 1.85 Gyr. This indicates 
that the lifetime of the spiral arm is very short (~ 120 Myr). 



The angular pattern speed measured from the snapshots 
in Fig. m (middle row) is plotted in Fig. \5\ It is seen to de- 
crease with radius, such that it almost equals the circular ve- 
locity of stars (Fig. [1] also shows circular velocity in angular 
terms at t = 1.77 Gyr for reference) in the disc at all radii. 
This indicates that this spiral arm feature is co-rotating, 
and is therefore unlikely to be long-lived. This is confirmed 
from the snapshots of our simulation shown in the top and 
bottom rows of Fig. [H which show that the arm starts de- 
veloping around t=1.73 Gyr and is winding and disrupted 
around t=1.85 Gyr. Hence the lifetime of the arm is about 



120 Myr and the arm is indeed transient. Bifurcations and 
breaks in the spiral arm features are seen regularly, which 
occur as the arms wind up. Wada et al. (2011) also find that 
the spiral arm features in their simulations are co-rotating, 
winding and short-lived. 



As an additional test, and to assist intercomparison 
with previous studies, w e compute spectrog:ra ms for m = 2, 3 
and 4 modes, following iQuillen et alJ (|201ll ). At each time 
step, which spans a time period of 1.28 Gyr centred on 
t ^ 1.77 Gyr, and each radial bin, we calculate 
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Figure 5. The pattern speed calculated between t= 1.75 and 
1.77 Gyr (dot-dashed red) and between t= 1.77 and 1.79 Gyr 
(blue dashed). The circular velocity at t= 1.77 Gyr shown in Fig. 
[T]is also plotted (black solid). The pattern speeds agree well with 
each other, and exhibit a decreasing trend that traces the circular 
velocity very closely over the radial range shown. Similar trends 
are found at other time steps with different arms. 



Wc{r, t, m) — cos {mOi 



Wsir,t,m) 



E 



sin (mOi 



(4) 



(5) 



where 6i is the angle at the position of the particles within 
the radial bin. We calculate the amplitude of each mode as 
a function of radius and find that all modes show similar 
strength at t = 1.77 Gyr, hence there is no one mode that 
dominates over the others. We then compute the Fourier 
transform in the above period. 



rT2 



[Wc{r,t,m) ^ iWs{r,t,m)y'"'h{t)dt (6) 



where h{t) is the Manning function. We compute the power 
at each frequency ranging from zero to the Nyquist fre- 
quency, which is shown in Fig. (6] The spectrogram anal- 
ysis shows us the significance of the pattern speeds of the 
wave modes (which we call mode pattern speed), which can 
be different from what is shown in Fig. [5] (i.e. the pattern 
speed of the spiral arm feature), if multiple wave modes 
interfere with each other. All three modes, especially the 
m = 2 mode, show that there is some power in the mode 
pattern speed that overlaps the circular velocity at many 
radii. There are also several horizontal features which could 
mean that there are a number of wave modes with con- 
stant, but different pattern speeds (Sellwood (2011) and 
references therein), that span different radii, and may be 
constructively and destructively interfering with each other. 
For example, the m — 2 mode could be interpreted as 
having two pattern speeds (^ 30 and 45 km s~^ kpc~^), 
with a faster inner pattern and a slower outer pat tern, as 
suggested bv iQuillen et al.1 (|201lh : ISellwood LinI ([l989h : 
iMasset Tagger! (|l997l l. Fig. [6] also indicates the 1:2, 1:3 
and 1:4 Lindblad resonances, i.e. = Q±K,/m, for m = 2, 3 
and 4 respectively. 



We note however some caution with regard to the spec- 
trogram analysis method when applied to transient, variable 
amplitude wave modes. We construct a toy model in which 
we set a base density of stars and distribute them randomly 
in a ring. We then add a small fraction (2.5%) of those in 
the ring at ^ = 0, tt to mimic the m — 2 modes. We impose 
a single, constant rotational frequency on the particles, and 
perform the same spectrogram analysis as above. Fig. [7] (top 
left panel) shows the case assuming 31 km s~^ kpc~^ rota- 
tion, for which this frequency is correctly extracted. The top 
right panel shows the case for two overlapping wave modes 
of 31 and 62 km s~^ kpc~^ rotation speeds, and both are 
shown to be extracted. 

However, the amplitude of the w ave mode can ev olve 
and/or disappear as demonstraced in ISellwoodl (|201lh . To 
demonstrate the effect of variable amplitude, we return to 
a single pattern speed of 31 km s~^ kpc~^, and repeatedly 
increase and decrease the number of m = 2 mode particles 
using a sinusoidal curve from to tt, that spans over 120 
Myr, which is roughly consistent with the lifetime of the 
arm feature in Fig. [5] The effect is shown (Fig. [71 bottom left 
panel) to create three peaks, with one at the input frequency, 
and two others: one either side of the real one. 

The appearance of new mode can shift the azimuthal 
position, which we model as randomly shifting the azimuthal 
position of the mode particles within to tt, every 120 Myr. 
This effect plus the variable amplitude is shown in the bot- 
tom right panel. This produces several peaks, none of which 
corresponds to the real input frequency. Although this is a 
simple exercise, and our aim is not to explore many possible 
cases of transient wave modes, this demonstrates that espe- 
cially for transient waves modes, there is a possible danger 
that the spectrogram may not show the real pattern speed of 
the wave modes. Moreover, any systematic variability could 
be responsible for the horizontal features seen in the spec- 
trogram. 

In any case, if there are indeed several wave modes 
present, it is evident that they must conspire in a specific 
way in order to produce a spiral arm feature that is 'appar- 
ently co-rotating' as is clearly shown in Fig. [5] Associating 
the multiple wave modes to the apparent spiral arm fea- 
tures is beyond the scope of this paper. We rather focus on 
how this apparently co-rotating spiral arm feature affects 
the stellar motion. Therefore, we discuss the remainder of 
the paper in terms of what we see in Fig. (5] We term this 
'apparent pattern speed' and 'apparent co-rotation'. 

A direct implication of a decreasing apparent pattern 
speed is that there are co-rotation radii all ove r the disc, and 
so th e radial migration (as described bv ISellwood BinnevI 
I2OO2I ) mentioned in Section 1 is expected to occur at a wide 
range of radii of the spiral arm. As described in Section 1, 
radial migration at the co-rotation radius has been predicted 
to preserve the circular motion of orbits (not to heat them 
kinematically) . 

In order to see how stellar motions are affected by this 
spiral arm whose apparent pattern speed deceases with ra- 
dius, we trace the motions of individual particles in our 
simulation, as in Fig. El First, we select a sample of par- 
ticles around our chosen arm of interest at 5.5 kpc radius at 
t — 1.75 Gyr when the arm is most prominent. The sample 
is selected to be near the plane of the disc (|z| < 200pc), 
and has a radial thickness of 0.25 kpc and azimuthal width 
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Figure 6. Spectrograms of the m = 2,3 and 4 Fourier components in a time window spanning 1.28 Gyr centred on f ~ 1.77 Gyr. 
The plots show the power in frequencies ranging from to the Nyquist frequency as a function of radius. Left: The spectrogram of the 
m = 2 component. Also marked is the circular velocity line (solid red), the inner and outer Lindblad resonances (solid blue) and the 4:1 
resonances (dotted cyan). Middle: The same as left panel but for m = 3 Fourier component, and 3:1 resonances shown (dashed cyan). 
Right: The same as left panel but for m = 4 Fourier component. 




Figure 7. Input pattern speeds of the m = 2 wave modes (dashed red line) and the power recovered from the spectrogram analysis 
(solid black line). Top left: Single Qp = 31 km s"-*^ kpc"-*^ is input and correctly extracted. Top right: Two wave modes of Qp = 31 and 
62 km s"-*^ kpc"-*^ are input. Bottom left: Same as top left, but the amplitude of the wave mode is varied. Bottom right: Same as bottom 
left, but random offsets of phase angle are implemented. See text for details. 
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of '^l radian centred on highest density point of the arm 
(left-middle panel of Fig. [8|). We follow their motion with 
respect to the spiral arm around which they were selected 
m Fig. El The middle row shows the snapshot and smoothed 
normalised density plots at the time they were selected. The 
preceding and succeeding rows show that of GOMyrs before 
and after respectively, and we move in a non-inertial frame of 
40kms~^kpc~^ from the middle step, hence only the middle 
row shows explicitly the actual co-ordinates in the smoothed 
plots. 

It is seen that as the spiral arm grows stronger, par- 
ticles from both sides of the arm begin to join the arm, 
which indicates the apparent co-rotation with the star par- 
ticles. We find that the spiral arm develops in a way 
akin to swin^ amplifica t ion (iGoldreich Lvnden-Bellll 19651 : 
ljulian Toomre(ll966l : iToomrd llQSli T In swing amplifica- 
tion theory, a leading arm can grow in density as it shears 
into a trailing arm. The star's epicycle phase and shearing 
motion of the arm conspire in a way that the stars in the 
spiral arm remain in the over-dense region for longer. This 
means that as the leading arm turns into a trailing arm ow- 
ing to the shear motion, self- gravity becomes stronger and 
accumulates more stars to the arm, and the amplitude of 
the arm begins to grow non-linearly. 

The left column in Fig. [8] indicates that the chosen par- 
ticles look like a leading feature at t=1.69 Gyr, which be- 
come part of a trailing arm at t = 1.75 Gyr as it wound 
up. The selected particles around the arm at 5.5 kpc (white 
dots) join the arm from the outer (inner) radii at the leading 
(trailing) side of the spiral arm while at the same time the 
spiral arm appears to grow in density. We therefore witness 
swing amplification in action. Note however that this process 
is different from the classic swing amplification mechanism 
where the co-rotation radius is assumed to be one specific 
radius. What we have found may be described as swing am- 
plification occurring over a wide range of radii. This is ac- 
companied by strong radial migration, which we describe 
below. 

From the selected sample of particles, we compute the 
angular momentum change, AL, over a period of 80 Myrs 
and choose those that exhibit the largest values of AL. As 
a fraction of their initial angular momentum, L, these have 
typical values AL/L 10 — 20%. We term these particles 
extreme migrators. The middle column of Fig. [8] shows the 
evolution of the extreme positive (cyan) and negative (pink) 
migrators. The 'positive' migrators are the particles that mi- 
grate towards the outer radii on the trailing side of the spiral 
arm. They are trapped by the potential of the spiral arm, 
which accelerates them. The co-rotating nature of the spiral 
arm feature ensures that during migration to outer radii, in- 
stead of passing through the spiral arm they remain on the 
trailing side (middle and right hand panels of Fig. [8]). There- 
fore, they continue to accelerate until the spiral arm is dis- 
rupted. The 'negative' migrators are particles that migrate 
towards the inner radii on the leading side of the spiral arm. 
They are decelerated as they become caught in the potential 
of the spiral arm, and because of the apparent co-rotation, 
they continue to decelerate on the leading side of the spiral 
arm until the spiral arm is disrupted. This illustrates the 
different motion that occurs on each side of the spiral arm. 

To demonstrate that these stellar motions and strong 
migration occur over a wide range of radii. Fig. [9] shows the 



same dynamical evolution as Fig. [8] for a sample selected at 
the radius of 7.5 kpc. This 7.5 kpc sample and the extreme 
migrators were selected using the same criteria as the 5.5 
kpc sample, however at t = 1.77 Gyr, since we find that 
the spiral arm at the outer radii grows later. It is clearly 
demonstrated here that exactly the same type of motion 
expected at the co-rotation radius happens at 7.5 kpc. This 
motion is at least consistent with the apparent co-rotation 
found in Fig. (5] because the particles continue to join the 
arm from both sides at a large range of radii as they migrate, 
instead of passing or being passed by the spiral arm, which 
is expected if the pattern speed is constant as predicted by 
density waves. This strengthens our conclusion of apparent 
co-rotation of the spiral arm feature made from Fig. \5\ 



3.2 Energy and Angular Momentum Evolution 

When discussing the changes in angular momentum associ- 
ated with radial migration, it is important to look at how the 
orbital energy of a star is affected. The relevant question be- 
comes: 'Is the star scattered during the migration by gaining 
significant energy of random motion?'. To shed light on such 
a question, we calculate the energy, and angular momen- 
tum, L, of the extreme migrators in Figs. [8] and [9] at 40 Myr 
before and after the time step at which they were selected, 
and show this in Fig. [TOl We call these two time steps the 
'initial' and 'final' time steps respectively. The solid black 
line indicates the L and E expected for a pure circular orbit 
at each radius. This represents the minimum energy which a 
star particle can have at a given angular momentum. In Fig. 
llOlwe show the position of the extreme migrators at the ini- 
tial (yellow diamonds) and final (red triangles) time steps. 
We can see that the negative migrators (top panel) and pos- 
itive migrators (bottom panel) move along the circular ve- 
locity curve in opposite directions to each other. Because 
they keep close to the circular velocity curve after migra- 
tion, their orbits must retain a very similar eccentricity, and 
so they gain little random energy and are not scattered into 
higher energy orbits (Sellwood &; Binney 2002). 

However there is some movement away from the circular 
velocity curve, and this corresponds to heating effects. We 
therefore test the degree of orbital heating of migrators with 
a larger sample of stars at different radii. 

We select the most extreme migrators from a sample of 
particles around the spiral arm in a radial range of 4-9 kpc 
at t = 1.77 Gyr. We plot in Fig. [11] their change in angu- 
lar momentum between initial and final time steps spanning 
a period of 80 Myr as a function of their angular momen- 
tum at the initial time step of t = 1.73 Gyr. We selected 
the particles with |AL| > 0.6, and calculate AE / AL for 
each particle, which we show plotted against the mean ra- 
dius of each particle over the 80 Myr period in Fig. [121 The 
filled (empty) circles represent the positive (negative) migra- 
tors. A clear decreasing trend with radius is seen for both 
sets of migrators. This trend qualitatively matches what is 
expected for AE / AL along the curve of circular orbit in 
Fig. IIOI but the overall position of the particles is slightly 
lower. This offset is evidence of the slight heating of the 
negative migrators. We note also that it is evidence for the 
slight cooling for the positive migrators. This is evident in 
the shallower gradients of particle movement in comparison 
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Figure 8. Left and Middle columns: the smoothed normalised density distribution (colour map) in the azimuthal angle - radius plane. 
Right column: the corresponding snapshots of the disc (the time referring to each row is indicated under each panel). In the left column, 
the particles (white dots) selected at 5.5 kpc at t= 1.75 Gyr are highlighted. In the middle and right columns the extreme migrators 
(see text) in the sample are highlighted by cyan (particles that migrate toward the outer radii) and pink (particles that move toward 
the inner radii) stars. Note that the coordinate of the angle in the left and middle panels at the top and bottom rows are shifted by the 
amount corresponding with Q = 40kms~^kpc~^, to keep the highlighted particles around the central region of each panel. 



with the tangents to the circular velocity curve for both sets 
of migrators. 

However the overall decreasing trend in parallel with the 
circular velocity line indicates that similar and only small 
heating losses occur over a range of radii. This demonstrates 
that radial migration can occur everywhere in the disc, and 
stars can slide along the spiral arms to different radii all over 
the disc with little contribution to disc heating. 



4 CONCLUSIONS 

We have performed N-body simulations of a pure stellar 
disc, and then performed a dynamical analysis of the spiral 
arms and particles around the spiral arms and traced their 
evolution. We come to the following conclusions: 



(i) We find in our simulations that spiral arms are tran- 
sient recurring features: we observe the continuous disap- 
pearance of spiral arms and the reappearance of new ones. 
This transient nature has always been found in numerical 
simulations, and is not consistent with spiral density wave 
theory. 

(ii) We performed two analyses on the pattern speed of 
the spiral arms: the stellar density trace method; and spec- 
trogram analysis. The latter shows a possible configuration 
of two or three wave modes that could mean that wave 
modes interfere constructively and destructively with each 
other. 

Regardless of whether this is the case, the stellar density 
trace method gave us clear evidence of an apparent spiral 
arm pattern speed that co-rotates with the stars at a large 
range of radii, which is reflected in points iii), iv) and v) 
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Figure 9. The same as Figure [S] but for a sample around 7.5 kpc radius. The particles were selected at 1.77 Gyr; one time step after the 
5.5kpc sample, in order we see the strongest, most clear radial migration. This is because the arm at smaller radii shows more prominence 
at earlier times than at larger radii. 



below . This is consistent with what is found in 'Wad a et alj 
(|201l h. and contrary to classic spiral density wave theory, 
which allows a constant pattern speed as a function of ra- 
dius, and hence there is only one co-rotation radius. Al- 
though we cannot reject that multiple modes are creating 
this apparent co-rotation, there is still the question of why 
the apparent pattern speed is always similar to co-rotation. 

(iii) Particles are shown to join the spiral arm from both 
sides at all radii. This is further evidence for the co-rotating 
nature of the spiral arm feature, because the arm must move 
at a similar speed to the particles in order for them to join 
the arm from both sides at all radii. 

(iv) Particles migrate radially along the spiral arm at all 
radii. Stars behind the arm are accelerated by the arm and 
slide along the arm to larger radii. Stars in front of the arm 
are dragged back by the potential, and slide down the arm 
to smaller radii. 

(v) Migrating particles do not actually cross the spiral 
arm. The co-rotating nature means that the particles stay 



on their respective sides of the spiral arm, so they are accel- 
erated (decelerated) until the spiral arm is disrupted. This 
means tha t radial migration is more efficient than that dis- 
cussed by ISellwood BinnevI (|2002h . with little contribu- 
tion to disc heating. 

Conclusions (i) and (ii) are already known from previous 
studies, and we include them to show the consistency of this 
study. Conclusions (iii), (iv) and (v) are new results which 
to our knowledge have not been reported before. 

A possible limitation of this study is that it focuses 
only on the stellar component. However, our preliminary 
work with simulations including gas, star formation and su- 
pernova feedback also show a consistent picture of transient 
co-rotating spiral arms and strong radial migration, which 
will be presented in a forthcoming paper. 

In this study, we have not addressed the mechanism 
of formation and destruction of the spiral arm structure. 
In this simulation, spiral arms appear to grow from small 
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Figure 10. The energy, E, and angular momentum, L, distribu- 
tion of the extreme migrators in Figs. [8] and [9] at 40 Myr before 
(yellow diamonds) and 40 Myr after (red triangles) the time step 
at which they were selected. The top (bottom) panel shows the 
results of the migrators that moved toward the inner (outer) radii. 
The solid black line indicates the circular orbit. Units are arbi- 
trary. 



shot (Poisson) noise in the (initial) density distribution and 
increase in prominence from there, possibly through mech- 
anisms similar to swing amplification. 

Although co-rotating spiral arms are expected to be 
short-lived due to winding, we observe that the spiral arms 
are disrupted before they are wound up very tightly. We also 
observe kinks i.e. change in pitch angle, in the spiral arms. 
The mechanism that controls the transient shape of the spi- 
ral arms and drives the disruption is unclear, and requires 
further work. 

In this study, we focused on the evolution between the 
eras of formation and destruction, and found that the spiral 
arms in our simulations are not consistent with the long- 
lived, rigidly rotating spiral arms of spiral density wave 
theory. The transient co-rotating spiral arm found in N- 
body simulati ons should be t ested by observations of both 
our Galaxy (ISellwoodI l201l | l and e xternal galaxies (e.g. 
Meidt et al.ll2009l : iFovle et al.ll201ll : ISpeights Westpfahll 
201lh . Although this paper focuses on the non-barred spi- 
ral galaxies, it is worthwhile to explore barred-spiral galax- 
ies where the spiral arms are affected by the bar (e.g. 




Figure 11. The angular momentum change over 80Myr for par- 
ticles in a radial range of 4.0 to 9.0 kpc around the spiral arm 
at t = 1.77 Gyr, as a function of their initial angular momen- 
tum at t = 1.73 Gyr. We define the extreme migrators to have 
|AL| > 0.6. The radius expected from circular motion with cor- 
responding angular momentum is also shown. 
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Figure 12. AE/AL for the extreme migrators selected in Fig. 
111! (open and filled circles) as a function of their mean radius 
between t = 1.73 to t = 1.81 Gyr. The filled circles are those 
with positive migration, and the empty are those with negative 
migration. The solid black line is the circular velocity. The units 
of y-axis or adjusted to match with that of Fig. [5] 



'Sparke k Sellwood" "1987"), and thus could be longer lived 
(e.g. Donner & Thomasson 1994 Binnev & Tremaine 2008; 
Baba et al.ll2009l . rQuillen et aLH201lh . We will study further 
the spiral arms in our simulations and endeavour to find 
observational consequences of the co-rotating spiral arm. 
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